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[Article by Katsuhisa Koura and Mikinari Takahira, of 
the National Aerospace Laboratory and Daiko Limited, 
respectively] 

[Text] Abstract: The “rarefied gas numerical wind tunnel 
(RGNWT)” constructed with the universal code 
described using the null-collision direct-simulation 
Monte Carlo method is used for the simulation of 
rarefied gas flows around the orbital re-entry experiment 
vehicle (OREX). Some aerodynamic characteristics and 
flowfield properties are presented. 


(RGNWT) which is designed for aerody- 
—FF oe 6 biol Gad Sah bo 
the space plane flies through rarefied gas. High- 
precision, high-efficiency simulation techniques are a 
must for three-dimensional full-airfoil simulations in 
rarefied gas. For this purpose, the null-collision (NC) 
direct-simulation Monte Carlo (DS-MC) method' is 

used where the mean molecular lower limit per collision 
cell is no more than approximately | and the computa- 
tion time is short. General-purpose use is made of a 
high-efficiency (uniform) rectangular parallelepiped for 
the computation area, collision cell shape, and data cell 
shape. The boundary conditions at the boundary of the 
computation area are that the uniform flow at the object 
upper flow and vertical and lateral boundaries are disre- 
garded, and the changes in the direction of flow at the 
object lower boundary are disregarded in making 


object are divided into shapes that are suited to the 
object shape independently of the collision cell shapes. 
As to the gas molecule reflection conditions at the object 
surface, a mixed form is assumed which is a mixture of 
mirror-surface reflection and diffuse reflection. 


The gas molecule model used is the variable soft sphere 
(VSS) model which was recently proposed as a high- 

precision, high-efficiency molecular model.” We are also 
now developing a molecular model for conducting high- 

precision, high-efficiency simulations of actual gas 
effects such as molecular spin, vibration, dissociation, 
electronic excitation, ionization, and radiation. 


2. OREX 


We used the rarefied gas numerical wind tunnel to 
conduct three-dimensional rare flow simulations around 

















Figure 1. Collision (Data) Cell and Object Surface Cell 





the orbital re-entry experiment vehicle (OREX). The 
collision (data) cells and the object surface cells are 
diagrammed in Figure l. 

Because the inter-molecular potential at high energy is 
indeterminate,’ we used values for the parameters in the 
VSS molecular model which correspond to those of rigid 
sphere models. We also assumed a uniform flow temper- 
ature dispersion wall for the object surface. The figures 
which follow depict the results with angle of incidence a 
= 0° for cases where flight is at a uniform flow Mach 
number M. = 27 (speed ratio S, = 25) in a rarefied flow 
with a uniform flow Knudsen number Kn,,, = 10, 1, 0.1, 
and 0.01, with the OREX diameter as “the standard 
criterion. 


The object axis (rotational) symmetrical surface infra 
density and temperature distributions are depicted in 
Figures 2 through 5 for Kn. = 10, 1, 0.1, and 0.01. Along 
with the reduction (high-level reduction) of the Knudsen 
number, the disturbance area in front of the object is also 
reduced. The density reaches several tens of times the 
uniform flow value in the vicinity of the forward stagna- 
tion point, while the temperature exceeds a hundred 
times the uniform flow value in front of the forward 
stagnation point. There is also a high temperature region 
in the back flow of the object which is several tens of 
times the uniform flow temperature. 


In Figure 6 are plotted the changes based on the Knudsen 
number Kn. of the drag coefficient C,. The drag 
coefficient declines along with the decline in Kn... The 
degree of convergence in the computational results for 
the decrease in the collision cell size is indicated on the 
diagram also, but the convergence is inadequate at Kn, 
< 0.1, and we see that it has become a largish value. 


The distributions of both the non-dimensional (surface) 
pressure P,/(p..¢ 


-’) on the OREX surface and the 
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Figure 4(a). Density Distribution; Kn... = 0.1 Figure 5(a). Density Distribution; Kn... = 0.01 























Figure 4(b). Temperature Distribution; Kn. = 0.1 Figure 5(b). Temperature Distribution; Kn. = 0.01 
































a 
2.5 age a 
Pas 4 | 
| | 
2.0 — —— 
Cy ey i 
: 
va 
0.01 0.1 J 10 & 
Ne 


Figure 6. Drag Coefficient; (Circle), Collision Cell = 
Computation Region/S0 x 50 x 50; (Square), Collision 
Cell = Computation Region/100 x 100 x 100 























Figure 7. Surface Pressure Distribution; (Square), Kn. 
= 10; (Circle), Kn., = 1; me Kn, = 0.1; + Kn, 
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The OREX simulations with the rarefied gas numerical 
wind tunnel were conducted in a joint research effort 
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[Article by Hajime Miyoshi of the National Aerospace 
Laboratory] 


[Text] Abstract: In this paper, it is shown that the 
Numerical Wind Tunnel (NWT) with an actual perfor- 
mance of more than 100 times higher than Fujitsu- 
VP400 is feasible using a multi-computer architecture 
and crossbar interconnection network. 


1. Introduction 


Let me tell you what I did at last year’s 3d aircraft 
computational aerodynamic symposium. | proposed: 


(R1) That it is necessary to quickly develop a high-speed 
computer which will exhibit at least 100 times the CFD 
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program speed of the Fujitsu-VP400 (herein- 
after called the VP400), with a main memory area 
of 32 GB or so for date storage, and to use this in the 
NAL computer system. 


This would be for the purpose of exploiting the merits of 
numerical simulation , for the purpose of 
employing the CFD gains made thus far in the aircraft 

t which is beginning to be done 
in Japan, and to push CFD to even higher levels. I also 


In this paper I press last year’s argument even further 
and report on: 


(1) the considerations that form our starting point for 
studying the idea of a numerical wind tunnel, 


ya aghe orp legge a Med egress com- 
puter exhibiting the seems and Li Oi) 
(3) the hardware for a computer exhibiting 


configuration 
the performance noted in (R1), and 


(4) a summary of the functions required in program 


Furthermore, with the purposes for developing a com- 
puter that will satisfy the performance demands noted in 
(R1) firmly in mind, in the remainder of this paper l 
consider that computer to be a numerical wind tunnel 
(abbreviated NWT). 


2. NWT Concept Starting Point 


In this section | will briefly discuss some aspects of the 
starting point for the NWT concept. The details are 
discussed in the References.” 


2.1 NWT Approach 
(1) NWT Will Be Distributed Main Memory Type Par- 
allel Computer 


In terms of how the main memory is structured, we may 
divide parallel computers into those which employ 


common main and those which employ distrib- 
uted main memory (cf. Figure 1). We decided to make 
the NWT a distributed main parallel computer 


memory 

for three reasons, namely (1) the fact that there are few 
limitations on improving raw computer processing 
speeds, (2) the fact that there are great possibilities for 
improving performance in the future by incorporating 
advances in element technology, and (3) the fact that the 
big problem with distributed main memory type parallel 
computers—i.e. the problem of addressing main mem- 
ory—can be to some extent resolved, we now aa 
incorporating virtual common main memory space 


(2) Element Computer (abbreviated PE) Will Be Vector 
Computer Equipped With Vector Registers (abbreviated 
VR) 
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Figure 1. Distributed Main Memory Parallel Computer. 





We have learned from our experience to date with vector 


computers that VR-equipped computers having chaining 
functions are suitable for CFD 


the vectoring ratio for CFD programs is 99% 


Paralleling and Vectoring Ratios Reacher 99% and 
Higher in CFD Programs. 


Most of NAL’s CFD programs now attain a vectoring 
ratio of 99% or higher. If the degree of paralleling is 


(1) Theoretical Elements 


su ded canted cate tad 0 toe 8 a Reclennea 
axis. In Figure 3 are plotted the trends in theoretical 
element densities, with gate density per chip plotted on 

axis and the year on the hori ontal axis. 
Bascd on these two figures, and on recent theoretical 
element per-chip power consumption trends, we expect 
that the following three devices can all be used as NWT 


(i) 30,000 gates/chip: Switching delay time 50 ps (10°? 
seconds); Power consumption 40 W/chip ECL gate array 


(ii) 100,000 gates/chip: Switching delay time 350 ps 
(10°'? seconds); Power consumption 20 W/chip 
BiCMOS gate array 


(iii) 200,000 gates/chip: Switching delay time 600 ps 
(10°'? seconds); Power consumption 10 W/chip CMOS 
gate array 


(2) Main Memory Elements 
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Figure 3. Theoretical Element Density Trends 





As to the main memory element, the 4 M/chip SRAM is 
beginning to be mass-produced. Accordingly, it is pos- 
sible to use either a 4 M/chip or 1 M/chip SRAM as the 
main memory elemen‘ and implement a main memory 
of 30 to 40 GB.' Whici\ chip to go with will depend on 
cost and speed considerations. 


3. Conditions Needed for in:olementing Numerical 
Wind Tunnel 


In the introduction we noted that the performance goals 
for the numerical wind tunnel have been established, but 
the answers forthcoming from computer engineering do 


namic testing equipment that complements the conven- 
tional wind tunnel. If the cost is too high, therefore, the 
NWT cannot be used in parametric studies and the like 


appropriate improving reliability and 
maintainability? The users of the numerical wind tunnel 
will be CFD researchers and technicians, rather than 
computer engineering specialists, so a stringent attitude 
will probably be taken towards malfunctions. Further- 
more, since the numerical wind tunnel is not just an 
experimental computer but also a practical computer, it 
will be just as important—in terms of the success or 
failure of the numerical wind tunnel—to maintain NWT 
hardware reliability at high levels as to achieve the 
targeted performance in the first place. 


3.3 Performance Demands for Element Computer 


Two items in particular will be necessary in order to take 
the NWT, which is a distributed main memory parallel 
computer, and achieve its widespread use among ordi- 
ony operators who are not computer specialists, as 
ollows. 


(i) Either an optimizing compiler for Fortran programs 
written according to a Fortran standard for parallel 
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programming, or a compiler that automatically opti- 
mizes programs written in current Fortran for parallel 
impiemcntation. 

(ii) Parallel-running programming support software 
including debugging and tuning tools. 


We will bring a practical computer in among CFD 
technicians and researchers who are not computer spe- 
cialists and seek 10 devise both (i) and (ii) above in 


conversations with and ordinary 
users. So doing, we will adopt a policy of promoting the 
growth and of knowledge needed for high- 


— coe 


(2) Research is being done on non-struciure lattices as 
well as on structure lattices. 


o> ae ee eee, ae ae 
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computation for single-axis directions in three- 
dimensional space. 


With structure lattices which are graphically simplified, 
it is not practicable to disregard the computation 


methods that employ sequential computations in single- 
axis directions, so, when one axis is used in vectoring, 


the remaining axis will be parallelized. 


When we consider that the scale of problems visualized 
for the NWT having the power set forth in (R1) is on the 
order of several thousand x several hundred x several 
hundred, this comes bouncing back in the form of a 
limitation on the number of parallel computers. 


In particular, with the IAF technique that is an impor- 
Ah A RAAT Dh ng 
t performance demands made on intercon- 
cedllia wines Gin wantin teh climes Gomente 
are also made on the network topology and on the 
number of element ‘omputers configuring the numerical 
wind tunnel. 
Our numerical wind tunnel! must be able to deal effec- 
tively with all CFD techniques. And we should not adopt 
an interconnecting network which exhibits high perfor- 
mance only for specific CFD techniques and inadequate 
performance for other techniques. 
From (R2) set forth in 3.3 above and the conditions 
noted in this subsection, we can see that it will be very 
difficult to implement NWT with several thousands to 
several tens of thousands of microprocessors. 


4. NWT Hardwar: Configuration, Performance 


4.1 PE Configuration, Performance 


The processing speed of a VR-equipped vector computer 
is determined mainly by the following eight factors. 

(i) Machine clock cycle time (hereinafter abbreviated 1) 
(ii) Pipeli tiplici 


(iii) Type and number of pipelines, and number of 
pipelines supporting simultaneous activity 


(iv) Ability of main memory to provide data 
(v) VR number and length 

(vi) Rise time for various pipelines 

(vii) Instruction execution control format 


(viii) Scalar instruction processing mode, register config- 
uration, cache memory volume and control mode. 


First, we assumed that items (vii) and (viii) would be the 
same as in the VP400, and then prepared many PE 
candidate models in terms of PE design parameters for 
the hardware resource quantities and timing noted in 
items (ii) through (vi). Next we examined typical CFD 
programs at NAL and selected 18 representative charac- 
teristics including number of vector instructions, ratio of 
various types of vector instruction, and the ratio between 
vector load/store instructions and vector operation 
instructions and then estimated 18 DO-loop processing 
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simulates the actions of a vector . The error in 
the VTAP simulations for the VP400 speeds 
for the 18 DO loops was within +/- 10% and the mean 
error was 6.8%. 

In Figure 4 we have plotted the results for one PE model 
series. The multiplicity in this PE model series 
was varied from | to 2 to 4 to 8 to 16, with the other 


parameters being made up of six fixed models. 


Pipeline multiplicity is plotted on the horizontal axis in 
this figure, and machine clock time t on the vertical axis. 


Here we modify the (R2) demand that the PE CFD 
program processing speed be VP400 or better to (R2’). 


(R2°) VTAP estimated value for PE speed for 
all 18 DO loops at vector lengths of 32 and 128 > VTAP 


lations on various other PE models (besides this P 
model series), we found that this PE model series is the 
most e<fective series among the PE models which we 


After comparing the power consumption of a PE with t 
= 3.1 ns and a pipeline multiplicity of 2 using ECL 
elements with that of a PE with t = 12.6 ns and pipeline 
multiplicity of 8 using BiCMOS elements, we found that 
the latter PE model would be better for use in the NWT. 





























When BiCMOS elements are used, moreover, one can be 
sure that a machine clock time of around 10 ns will be 


In Table | [not 


| 
| 
: 
i 


An order to satisfy R1, it is only necessary to have 150 
or more PE units. If we have 256 MB main memory 
capacity for each PE unit, then the total main memory 
capacity for the NWT would be about 38 GB. 


In Figure 5 is presented a conceptual block diagram of 


5. Features Needed in NWT Programming and 
Compiler (P&C) 


In this section we give a broad overview of the main 
P&C features required for maintaining an NWT PE 
efficiency of 0.5 or better. These features may be imple- 
mented in program writing (the programmer's responsi- 
bility) or by the compiler, or by a combination of both. 


In Figures 6-1 and 6-2 are given conceptual running 
diagrams when programs are run on one PE and on n PE 


units, respectively. 
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Figure 6-1. Conceptual Diagram of Executing Job With 1 PE Unit 
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Figure 6-2. Conceptual Diagram of Executing Job With n PE Units 





S,: Execution time required for i’th sequential process in 
program 

SO,: OS execution time required for S, p:ocess 

SI,: OS execution time for I/O processing required for S, 
process 

I/O,: VO Execution time required for S, process 

P;: Execution time when j’th parallel-processable portion 
in program processed by | PE 

PO;,: OS execution time required for P, process 


PI,: OS execution time for I/O processing required for P, 
process 


I/O,: ./O execution time required for P, process 
ParaOS,: OS processing time that becomes necessary to 
do parallel processing 


P,,: When 


P,, execution time for P, 
potion tat shouldbe proceed by PE; (2 by the 1st 


S/RO,,: OS execution time required to S/R (data 
transfer) data needed for P,, process 


S/R,,: Time required to send data 


PO,,: Execution time for OS processing required for P, 
process 


PI,,: OS execution time for I/O processing required for 
P,, process 


1/O,,: I/O execution time necessary to execute P,, 


sy Execution for OS processing required to synchro- 


, with other PEs 
W,,: PE, synchronization wait time = Execution time 


é. - SO; US; USI, UIA, 


P; 


when P, is parallel-processed 

The time (TS) for processing the program with one PE 
unit and the time (TPN) for ing the with 
n PE units are given by the following formulas. 


" 
TS = ~ §, + 3 {P, + PO; +P, +1, ) 
i we . 
TPN = ~ {S, + ParaOS, }+ ~ P, 


—AAMAM 2 ee oe 
operations (included only in P, or S,) required to execute 
the entire program, we get F/TS as the program pro- 
ents ek One Ot olen Cee 

the program processing speed with n PE units. Accord- 
ingly, it is essential to make TS/TPN as large as possible 
in the interest of improving speed by paralle! processing. 
The P&C features needed, in conceptual terms, to 
accomplish are evident from Figures 6-1 and 6-2. The 
most important ones are noted below. 


(A) Improved ability of P&C to find those portions that 
are parallel-processable so as to shorten S, and reduce the 


ParaOS, frequency 

















re SS. Seen Se eee 
data transmission volume and frequency so as to reduce 


frequency of S;, and shorten processing time 


(C) P&C functions to facilitate P;, and 3,, parallel 
action 

(D) P&C functions to facilitate Pl,, and P,, parallel 
actions 

(E) P&C functions ‘o facilitate reduction of volume and 
frequency of referencing data stored in other PEs for the 
purpose of reducing the SYO,, frequency 


(F) In the interest of decreasing W;, as much as possible, 


(i) P&C functions to facilitate equalization of divisions 
of P; to P;, 


(ii) Since it is necessary to make it so that S,, processing 
does not compete on the network more than the absolute 
minimum, in order to make S,, equal for all 1's, the 
necessary P&C scheduling functions 


There must naturally be cooperation by both the hard- 
ware and the computation methods to facilitate the 
achievement of (A) through (F). If the degree to which 
(A) through (F) are achieved rises, it will be possible to 
satisfy {R1) with our NWT CFD program processing 
speed using P&C. 
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Three Dimensional Calculation of Flow Inside 
Supersonic Inlet 
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[Article by Akira Fujimoto of Gifu Technical Institute, 
Kawasaki Heavy Industries, Ltd., and Keisuke Sawada, 
Institute of Fluid Science, Tohoku University] 


[Text] Abstract: A three-dimensional Navier-Stokes cal- 
culation of flow inside a supersonic mixed-compression 
inlet, designed for Mach 2.5, was carried out. A multi- 
block method was applied to the calculation grids. Com- 
plex configurations of bleed holes and throat slot were 
simulated in order to compare the solution with experi- 
mental (test) data. As a result, detailed information on 
the flow—which is usually difficult to obtain experimen- 
tally—was obtained. 
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At this inlet under the supersonic flow-through condi- 
the flow was found to be two-dimensional in 90% 
wise direction. In the vicinity of the side- 


SS ee oe S ee & 
applications in such vehicles as the next generation of 
airliners and space planes.’ Research to date 


sonic mixed-compression air inlets using experimental 
techniques’ and two-dimensional mathematical 
analysis.’ The experimental techniques do provide much 
useful information, but there are limits as to model size, 


this information in subsequent research 

2. Calculation Grids 

For the calculation grids we used multi-block mapring. 
With the multi-block method, grids can be generated 
relatively easily for even the most complex shapes, and 
the analysis codes which use this require no special 
handling for each object shape. In Figure | is shown the 
relationship between objects in the grid space defined 
regularly by I, J, and K, and objects in actual physical 

(x, y, z). In Figure 2 are 


ii 
7 
is 
a8 


has a height:width proportion of 2:3, with the outlet 
shape not choked off in a circular shape but made a 
square shape (and the same flow path shape as was used 


in experiments described subsequently). 
3. Analysis Code 


The basic equations used are the Reynolds mean Navier- 
Stokes integration equations. Coakly’s q- model is used 

















(a) Overall View 


(b) Enlarged View 


Figure 2. Calculation grids 
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for the turbulence model.‘ Separation is determined by 
the limited volume method and cell boundary numerical 
flux by Roe’s approximating Rieman analysis, to which 
TVD conditions are applied.* For this calculation, the 
ADI method is used for the y-z plane and the Gauss- 
Seidel method for the x dimension. The equations are 
spatially solved, with the first-order-precision negative 
solution method used for the time integration. 


4, Analysis Conditions 


We made the Mach number 2.5, the Reynolds number 
10’ based on the standard of inlet height of the item 
being tested, with a = § « 0°. In this analysis, we made the 
area of the exhaust louver-type exhaust port constant 
without particularly regulating the bleed volume. In 
these tests, moreover, since the objective was to ascertain 
the interference in the supersonic section between the 
sidewall boundary layer(s) and the diagonal shock 
wave(s), we did not generate a terminal shock wave, but 
employed conditions such that outflow was still at super- 
sonic speeds at the main duct outlet. 


5. Analysis Results 


In Figure 3 [not reproduced] is represented an equal 
Mach contour diagram on the center surface. The two 
ramp shock waves and cowl shock waves, and the 
entropy compression continuing therefrom, were 
attained as designed, and the results obtained were 
almost the same as those from the two-dimensional 
analysis.’ The throat shock wave is a diagonal shock 
wave that is generated when the bleed causes the flow to 
strike the back of the slot. Since we have set low 
back-pressure conditions, the flow is accelerated behind 
this and flows out at supersonic speed. In Figure 4 [not 
reproduced] is shown an equal Mach contour diagram at 
the 95% span position. In this way, in the vicinity of the 
sidewall, we find separation in the vicinity of cowl shock 
wave injection point. The third-column ramp bleed 
hole(s) prevent this separation region from spreading, 
indicating the importance of bleeding at this position. As 
a result of varying the cross-section delineated by the 
equal Mach contour diagram in the span direction, we 
found that two-dimensionality is maintained to approx- 
imately the 90% span position. 


In Figure 5 [not reproduced] we have an equal Mach 
contour diagram at the x = constant cross-section as seen 
from different perspectives, (a) and (b). (The wall surface 
on the near side is not represented.) The separation 
region in the vicinity of the cowl shock wave injection 
point is sucked into the aft throat bleed hole(s), and we 
can observe a thick boundary layer region in the sidewall 
cowl region(s) also. 


In Figure 6 [not reproduced] we see flow trails from 
particles emitted from three places very close to the 
sidewall as seen from (a) the side, (b) diagonally from the 
front, and (c) directly from the front. Looking down- 
stream, the three-dimensional separation region at the 
left sidewall becoines a vortex which rotates clockwise 
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along the cowl shock wave and is eventually sucked into 
the throat bleed place. The flow which passes through the 
cowl shock wave is subjected to a sudden reverse- 

pressure gradient. The inner layer region of the sidewall 
boundary layer has lower kinetic energy than the outer 
layer, so the flow turns toward the ramp side (up). That 
is why a clockwise vortex is generated. 


When we view the movement of this vortex in moving 
pictures, we find that the speed of rotation increases the 
closer the vortex approaches to the throat bleed hole(s), 
with the occurrence of vortex stretching. 


In Figure 7 we compare the wallstatic pressures." At (a) 
is shown a static pressure distribution in the flow direc- 
tion, with made between the values on the 
center surface in this analysis (solid curve), the two- 
dimensional NS analysis values (dotted curve), and the 
test values (symbols).*:’ In the external compression 
region (x < 1.7), the two-dimensional and three- 
dimensional values agree perfectly, but a slight differ- 
ence develops in the internal compression region (x > 
1.7). We think this is due to the influence of the sidewall. 
Nevertheless, there is relatively good agreement between 
the two-dimensional and three-dimensional analysis and 
test values. At (b) in Figure 7 we have diagrammed the 
es ee i a ee ee te te 
positions indicated in (a) (i.e. Sec. 1, Sec. 2). The 
agreement with the test values is good. 


6. Concluding Remarks 


The following facts were established as a rent of con- 
ducting the three-dimensional analysis on the supersonic 
uir inlet. 

(1) With this air inlet, under conditions of * 


outflow, the two-dimensionality of the flow is main 
tained to approximately 90% in the span direction. 


(2) In the vicinity of the sidewalls, three-dimensional 
separation occurs along the cowl shock wave, and the 
third-column ramp bleed helps prevent this region from 
spreading. 


(3) With respect to the static pressure dis'ribution at the 
wall in the center plane under conditions of supersonic 
outflow, there is almost no difference with the [results of] 
two-dimensional analysis, and there is good agreement 
with the test values. There is also good agreement 
—“ the spanwise static enna oma and the 
test values. 


With these findings in hand, we plan to conduct three- 
dimensional analyses in cases where the flow is deceler- 
ated to subsonic speeds by the terminal shock wave. 
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Figure 7. Wall Static Pressure Distributions (Compared 
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[Text] Abstract: Numerical —*8 of .. 

averaged Navier-Stokes equations 

dimensional flows cae t aie > —— super- 

sonic air intake with its operation range from 
to subcritical. 


Computed wall pressure data are compared with wind 


operation is normally done in a supercritical state that is 





























near-critical, with a margin allowed prior to subcritical. 
When this is done, there is interference between the 
boundary layer(s) and the final shock wave inside the 
duct and this can trigger separation. If the vortices 
induced by such separation become large enough, not 
only is the flow pipe constricted so that the predicted 
pressure restoration cannot be effected, but the final 
shock wave may be pushed out so that a subcritical state 


PSs 


provided fur the intake outlet so that any desired back- 
pressure can be simulated. 


The main flow condition is sea-level Mach 2. 


4, 


The basic equations used are the two-dimensional Rey- 
nolds-averaged Navier-Stokes equations in which a two- 


layer algebraic turbulence model is employed. ' 


We perform spatial integration using the finite volume 
method, and implement second-order precision with 
MUSCL. For the evaluation of convective factor flux at 
the cell boundary we use Roe’s te Rieman 
solver? The factor is discretized with the 
second-order precision central differences. 


Time integration is performed with the second-order 
precision rational Runge-Kutta method.’ 


5. Boundary Covrditions 


As to the exit conditions inside the intake, supersonic 
i the virtual 











vortices, and inevitably a non-physics result 
is imparted. With the virtual nozzle the back pressure 
can be controlled and it is possible to set simple 
boundary conditions which make sense in terms of 


6. Comparing Calculation Results Against Test Data 


poe S at (a) through (c), are shown density distribu- 
at supercritical, critical, and subcritical states. As can 


strong, 

tices are large, while at critical the vortices do not grow all 
that The static pressure distribution at the lower 
wall is depicted in Figure 4. In this figure the test 
results for the supercritical state are represented by small 
circles. Except for the pressure rise caused by the second 
shock wave, the calculation results agree well with the test 
results. The calculation results coincide with the two- 
dimensional analysis solutions, so the differences in these 
results are thought to be due either to three-dimensionality 
or to measurement errors. The total pressure distribution 
inside the intake is di in Figure 5. It is evident 
that separation-vortex-induced distortion occurs close to 
the wall surface. 


We found that, in the Eulerian calculations that were 


method.* We, therefore, added other schemes and did 
the calculations using six different methods, and then 
compared the results. 


The first schemes we took up were the Roe method, van 
Leer method, Osher method,’ Harten-Lax-van Leer 
(HLLE) method,® and modified Harten-Lax-van Leer 
(HLLEM) method.® With the Roe method, furthermore, 
we employed the Chakravarthy’ and Lin® techniques in 
order to satisfy the entropy conditions. With the Roe 
method, it is necessary to add an artificial viscosity to 
the characteristic velocity A* in order to prevent expan- 
sion shock at points of sonic velocity. However, with the 
Chakravarthy technique, the A**”,,,,. at the cell 
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(a) Supercritical Mode 












































(c) Subcritical Mode 
Figure 3. Density Distributions 





obtained by Roe approximation is varied as 
follows, namely, at sonic rarefaction (A*, < 0 < A*,,;), 


Uep= tat Hak 


With the Lin technique, there is a converging accelera- 
tion effect, and an effect which inhibits the carbuncle 
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Figure 4. Static Pressare Distribution at Lower Wall 
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Figure 5. Total Pressure Distribution Inside Duct 





by which the Mach lines are distorted at 
the nose in hypersonic calculations for blunt objects. We 
next represent the Lin method. 
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Di, jor +2Di, j +Pi, j-1 
where p, , is the pressure. At the cell boundaries, we use 
averages of the values of adjacent cells. The values of the 


constants are k, = 0.25, 0 < k, <5, and k, = 15. Here we 
used k, = 0.25, k,= 5, and k, = 15. 


We show the results of Eulerian calculations done in 
relatively coarse meshes. The tip of the intake lip has 
become a flat head having twice the thickness of the grid 
interval in closest proximity to the wall (cf. Figure 6). By 
this means we obtained interesting results. In Figures 
7(a) through 7(f) are depicted the density distributions 
based on the Roe (Chakravarthy), Roe ** Osher, 
HLLE, HLLEM, and van Leer methods. The respective 
lower wall static pressure distributions are plotted in 
Figure 8. We can divide these into three ou 
according io the separation conditions, namely (1) Roe 
(Chakravarthy) and Osher, (2) Roe (Lin) and HLLEM, 
and (3) HLLE and van Leer. In group (1), separation was 
caused by the final shock wave, and the induced vortices 


pushed the shock waves upstream, resulting in subcrit- 
ical states. In group (2), slight separation does occur, but 


and %b). Inside this high-entropy layer, the velocity is 
slower than outside of it, resulting in a boundary layer- 
like velocity distribution. We think that separation 
occurs here due to a steep pressure gradient produced by 
the final shock wave. The generation of the high-entropy 
layer occurs when the tip angle is exceeded, and we do 
not think it occurs when the tip is made sharper. The 
group (2) high-entropy layers are thinner than those in 
group (1), for which reason we do not think that large 
separation was reached. In group (3), conversely, the 
high-entropy layer is thick and the peak is low. We 
believe this is due to the fact that the linear wave 
resolutions are inferior with scheme (3), as they are 
known to be. Thus we think that the differences in local 
linear wave resolution cause large changes in the numer- 
ical solutions. 


8. Concluding Remarks 


We were able to simulate flows in each mode, without 
going to extremes, by using a virtual variable nozzle to 
control back pressure. It is also possible to do unsteady 
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Figure 6. Intake Eulerian Calculation Grid 
(Lip Portion) 























Figure 7(a). Density Distribution Based on Roe 
(Chakravarthy) Method 























Figure 7(b). Density Distribution Based on Roe (Lin) 
Method 
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Figure 7(c). Density Distribution Based on Osher 
Method 
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Figure 7(f). Density Distribution Based on van Leer 
Method 
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Figure 7(d). Density Distribution Based on HLLE 
Method 
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Figure 7(e). Density Distribution Based on HLLEM 
Methoa 
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Figure 8. Lower Wall Static Pressure Distribution (Eal- 
erian Calculation) 
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clement which must handle either most or all of the 


putation technology, and the range of applications 
broadens year after year. In the field of compression 
flows, in particular, it has become possible to do high- 
precision, stable calculations for high-Mach flow sites, 
thanks to the development of the TVD scheme. As a 
result, it has become relatively easy to analyze the flows 
inside supersonic and scram-jet intakes, with an 
increasing number of such analyses being done with each 


passing year. 
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Romp dived Cavity bleed for the viscosity factor has been thin-layer approxi- 
__ Sedbsentic éiffeser mated. This formula can be expressed as foll | 
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Figare 1. Simplified Diagram of Sepersonic latake Test Q, & F 

ip) 
In this research, we analyzed the flows inside supersonic = ted | 
intakes having a cavity in the throat section, using the Q=v 
TVD scheme and solving numerically with the Navier- vy 
Stokes equations. We were thus able to simulate the flow ‘a4 
conditions in the vicinity of the cavity under various 
conditions. We also compared the results of these calcu- ee ae 
lations with the results of wind tunnel tests done on . 
models of the same intakes.' Bey": pul +£,p 

puU+&,p 
2. Basic Formulas, Numerical Calculation Methods (e+p)U | 

j ' 
2.1 Intake Model ev 

puV+ 
In Figure | is presented a simplified view of the super- Peay": wal | 
sonic intake used in the wind tunnel tests selected for our puV+n,p 
computation subjects. This intake is called a mixed- + — 
compression intake which uses both external and | (e+p) 
internal compression. It performs external compression , 0 ’ 
with four shock waves which are produced from the 
ramp, and then performs internal compression by three 2.2, Ou ws du av 
more shock waves from the cowl. The flow passes w (ni +05) =—+=05 (02 = +2 =) 
through the final shock wave of the throat, is decelerated ' — — o 
to subsonic speed, and is pressure-restored by a diffuser. §=s"; 2..2,0U A Ou du 
In a mixed-compression intake it is necessary to hold the a (nb +93) 5 + 519 1s 5 +995) 
shock wave(s) inside the channel. This means that we vo is 
may expect large separation to occur with it being a) #7 024,2)96 
Sanaee and tee Denier larer, P 2328 —— 

⸗ . For reason, . 

holes are positioned in the wall just in front of the 
—————— of the beondary layer. The where p is the pressure determined from state equations 
cavity that is made in the upper wall at the throat is for an ideal gas, Uand V are contravariant velocities, y is 
called a throat cavity and functions to hold the final the specific heat, Pr is the Prandtl number, Re is the 
shock wave stably in the throat. From this cavity there is Reynolds number, a is the speed of sound, and J is a 
a flow path which extends through to the outside of the | Conversion Jacobian. When doing the calculations, for- 
intake. Part of the flow inside the cavity passes through mula (1) is time-integrated with the cross-braced ADI 
this path and flows out of the intake. A flow plug that is method,” and steady solutions are obtained. The TVD 
provided at the diffuser exit moves back and to scheme of Harten-Yee° is used in discretizing the con- 
control the flow volume and thereby adjust the exhaust vection factors on the right side at this time, so as to 


2.2 Basic Formulas and Algorithms lence model.‘ The cavity and external flow regions are 


The basic formulas used are the two-dimensional com- _ calculated separately, but handled by expanding them 
pressible Navier-Stokes formula in which an ordinary _into a single calculation region. For this reason, we wind 
coordinate system (€, n) is employed. The formula used up having grid points in the solid regions of cowl and 
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Figure 2. Calculation grid (530 x 250) 





conditions of 2 x 10°. The conditions of the flow plug at 
the exit of the diffuser are represented as Ae/A 
area 


In Figure 3 are represented the results of computations at 
Ae/At = 1.13. At (a) the overall flow conditions in the 
intake are depicted in terms of a Mach n 

tion. At (b) we have an enlargement of 


head for the cowl tip. Then from the cowl we have three 
shock waves which extend as far as the throat cavity to 
produce the shock wave pattern as designed. The weak 
swell that is seen at the ramp wall just in front of the throat 
cavity is due to bleeding. We can see a slightly strong swell 
developing from above the cavity, and thus know that the 
flow is circulating inside the cavity. Part of this circulating 


i 
3 
§ 
2 


flow passes through the duct which goes to the outside of 
the intake and is thus bled off to the outside. 
compression wave and swell are vi 
the cavity. At (c) is seen a Schlieren photograph 

the vicinity of the throat in the tests. In this photograph 
can see the shock wave(s) coming from the cowl and 
way in which they reach as far as above the cavity, as 
as the compression wave and swell from the back end 


: 

| 

i 
ẽ 

E> 

Het ly 


beyond the end of the throat, so the optimal conditions are 
not attained. 


In Figure 4 are plotted the results of calculations at Ae/At 
= 0.94. First of all, from (a), when we look at the condition 
of the flow on the upstream side from the throat, there is 
no change from the computed results under the conditions 
previously indicated. This is because the entire flow 

from the throat is supersonic, and because the 
effects of the flow plug change are not propagated this far. 
Meanwhile, a terminal shock wave appears a little behind 
the throat cavity. This was not seen in the previous case(s). 
This shock wave appears to interfere with the compression 
wave from the back of the throat cavity and curve in an 
arc. In the Schlieren photograph at (c) we see a curved 
terminal shock wave a little toward the throat cavity. 
However, as may be seen from the computational results, 
this shock wave itself does not curve, but it appears to 
curve due to the interference between the terminal shock 
wave and the shock wave from behind the throat cavity. 
Compared to the previous case(s), the influence of the 
presence of this shock wave is great. Based on the wall 
pressure distributions plotted in (d), the pressure rise that 
occurs near the throat is far higher than in the previous 
case(s) when no terminal shock wave existed. 
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(b). Mach number contours around the throat. 
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(d). Static wall pressure distributions. 


Figure 3. Calculation Results (Ae/At = 1.13) 
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(b). Mach number contours around the throat. 
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(d). Static wall pressure distributions. 
Figure 5. Computation Results 











— So Sry hain = 0.86 where 


that this shock wave appears farther forward than in the 
F-evious case(s), and the shape thereof is not really 


ramp and cowl wall(s), even at the side walls. These 
eddies reduce the intake efficiency, and adequate pres- 
sure rise is not achieved. 


In Figure 6 is represented the relationship between 
captured flow volume (Gco) and total pressure recovery 
rate (Pr). The captured flow volume is roughly constant 
up to the point where the flow plug is closed to a certain 
point, and the total pressure recovery rate increases as 
the flow plug is closed. 


The captured flow volume decreases from certain condi- 
tions because the terminal shock wave reaches a point 
above the throat cavity, the flow downstream from the 
shock wave is pushed inside the cavity, and this flows out 
to the outside. The test results are also noted in this 
graph. When we compare these to the results, 
we see that there is qualitative agreement in the trends, 
but the computed values for both the total pressure 
recovery rate and captured flow volume are higher. We 
can say nothing for sure about the reason for this, but a 
disparity appears in the pressure distributions due to 
three-dimensionality, and the cavity bleed hole is three- 
dimensional. Hence the possibility that this also is due to 
three-dimensionality is great. 


4. Summary 


We simulated the flow inside a mixed compression type 
supersonic intake having a cavity in the throat section, 
using the TVD scheme to solve the Navier-Stokes equa- 
tions. In the analysis, we handled the intakes of complex 
shape in single computational regions. This enabled us to 
perform highly efficient calculations. In making compar- 
isons with test results, we found that the numerical 
calculations represent well the changes in flow condi- 
tions due to opening and closing the flow plug, thus 
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presence separation 
— Accordingly, if more precise intra-intake 
need to be predicted, it will probably be necessary 
to go to three-dimensional analysis. 
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[Text] Abstract: In this paper, th roles of wind tunnel 
testing and computational fluid dynamics (CFD) in the 
aircraft development is briefly described. The stress is 
put on the importance of CFD as an aerodynamic design 
tool to drastically reduce both cost and period of devel- 
opment. 


Prior to the application of CFD to practical design 
problem, the existing CFD codes must be verified with 
the reliable experimental data. A demonstrative example 
of CFD validation is presented to show the state of the 
arts of CFD by comparing the Navier-Stokes simu'ation 
results with the transonic wind tunnel test data for a 
NACA0012 airfoil at several angles of attack. 


The importance of application of CFD to inverse design 
problem is also discussed. 


1. Role of CFD in Aircraft Development 


It is no exaggeration to say that the wind-tunnel test used 
to be the only design tool available for aerodynamic 
design work. We create an artificial air flow inside the 
wind tunnel at the same velocity as during flight, we 
insert a small model of analogous shape to the actual 
craft, and we measure the pressures on the surface of the 
model and the other aerodynamic forces which operate 
on the model as a whole, thus obtaining the data neces- 
sary for the design. However, wind-tunnel tests have to 
be repeated time and time again until an airfoil shape is 
obtained which satisfies the desired performance figures. 
This requires the expenditure of enormous amounts of 
time and money. Furthermore, due to the interference 
caused by the wind-tunnel walls and model supports, and 
to the Reynolds number effect, it is very difficult to 
create the same conditions of air flow inside a wind 
tunnel as occur in actual flight. As was pointed out early 
on, the influence of the Reynolds number becomes 
increasingly problematic in aerodynamic design as air- 
craft become larger and faster. 

But things have changed dramatically in the field of 
aerodynamic design in recent years with the advent of 
the supercomputer and the amazing ts in 
CFD (computational fluid dynamics). By leaving to CFD 
many of the things that previously had to be done in 
wind-tunnel tests in the aerodynamic design process, the 
cost and time required for model fabrication and wind- 
tunnel operations can now be drastically reduced. CFD 
is particularly well suited to parametric studies involving 
shape variation and the like. CFD is expected to play an 


At the current time, however, due to the computing 
times involved, almost all of the computer code used in 
design work that is three-dimensional is potential flow 
analysis code.' As a consequence, it is very difficult to 
evaluate aerodynamic characteristics at off-design points 
where viscous effects are prominent. So wind-tunnel 
tests are also required in the preliminary design stage. In 


order for Navier-Stokes code to be used freely as the 
design code in actual design work settings, the entire 


operation—beginning with airframe shape definition, 
including grid generation and numerical computation of 
flows, and ending with the derivation of the final aero- 
dynamic coefficients—must be concluded within a single 
dat,” plus the computation time for each case for the 
same shape must be on the order of several minutes. 
Essential to achieving this will be the development of a 
numerical wind tunnel having functions for the auto- 
matic generation of airframe surface grids and space 
grids, and a computer with the power to process at 100 
GFLOPS or faster. 


At the stage of the main design or detailed design, 


2. Wind-Tunnel Tests and CFD Reliability 


In Figure 1 we have a single diagram that summarizes 
the results of testing the NACA0012 wing shape in a 
typical two-dimensional wind tunnel.ꝰ The shock wave 
position (distance from leading edge taking wing chord 
length as 1) is plotted on the vertical axis, an 


attack, and wing symmetry, the position of the shock 
wave is greatly different from wind tunnel to wind 
tunnel. A number of factors may cause these differences, 
including interference between the top and 
left and right wind tunnel walls, air flow bending, error 
in Mach number setting, model fabrication errors, and 
measurement errors. By looking only at this diagram 
there is no way in the world for us to determine which 
data are more correct, but McCrosky’ analyzed these 
wind-tunnel test data and predicted that the true shock 
wave position must be somewhere between the two 
dotted lines that are parallel to the horizontal! axis (i.e. 
between 44% and 48% of the chord length). The data of 
Harris (NASA)—which is considered the most reliable 
he omens data and is widely cited to demonstrate 

in on the upper dotted line (48% 
line, corrected [oquare) marks), while the data fro the 
NAL two-dimensional transonic wind tunnel ([bullet] 
marks, corrected for wind-tunnel wall interference)‘ fall 
in between the dotted lines, even though showing some 
variation. (Note that when no corrections are made 
({circle}, [triangle] marks), the shock wave is positioned 
considerably upstream.) 


As seen in the foregoing, it is possible to considerably 
reduce the variation between different wind tunnels by 
correcting for the wind tunnel walls, but it is utterly 
impossible to make all the data agree. Each wind tunnel 
has its own peculiarities, and it is next to impossible to 
eliminate these. Not only this, but when the same model 
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@ Harris, no trip 
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Hacris, corrected by Sewall 

Vidal, AR=8; trip 
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D Ohman, side-wall suction, no trip 
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@ Thibert, porous walls, AR=2.7; no trip 
@ Lee & Gregorek, porous walls, AR=1; no trip 
O Kraft, adaptive porous walls, AR=3; no trip 
@ Gregory & Wilby, slotted walls, AR=1.4; trip 
 Navier Stokes calculations, fully turbulent 
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Figure 1. Shock Wave Position vs. Reynolds Number 





is retested on a different day, there is usually some 
disparity between the two sets of data. The problem is 
so-called reproducibility. What then about CFD in this 
respect? The situation is absolutely the same. There will 
be as many answers as there are people making the 
computations, as there are schemes «sed, and as there 
are grids. The designer does not know which data can be 
relied upon. Well, at least this was the situation in the 
early days of CFD development. But CFD has developed 
rapidly since those early days, and its effectiveness is 
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now as widely recognized as its limitations. When there 
is either very little or no separation in the flow, so long as 
we use schemes and turbulence models that are reliable 
in the foregoing examples, for instance, almost all of the 
solutions are nearly the same if the number of grid points 
is large enough. And when it comes to reproducibility, 
CFD is perfect. 


We next turn our attention to a comparative study of 

generated numerical solutions and wind-tunnel 
data on the NACA0012 airfoil. The basic equations used 
in the computations were the thin-layer-approximation 


Navier-Stokes equations and, for the 
turbulence model, the Baldwin-Lomax algebraic model.° 
The scheme was the TVD difference method® 


grid points in the flow direction and in the direction 
perpendicular to the wing plane, respectively, were 

increased in —3 as soon as the drea coulficlens 
ceased to vary that was determined as the optimum grid 
(approximately 25,000 total grid points with this solver), 


In Figure | are noted both the computed shock wave 
positions for Re = 2.1 x 10’ and 1.2 x 10° ({diamond] 
marks). The shock waves are positioned at 48.5% for 
high Re numbers and at 47% for low Re numbers. From 
these results, we think that the true position of the shock 
wave is somewhat downstream of McCrosky’s predic- 
tion. We next compare the NAL two-dimensional wind- 
tunnel data, Harris data, and McDevitt data with the 
pressure distribution at Mach = 0.8 and Alpha 


believe this is either due to the wind-tunnel model being 
made incorrectly or a failure to completely eliminate the 
side-wall effects. In actuality, as is indicated in Figure 3, 
when the model fabrication precision was enhanced, and 
the tests were run again with the model’s aspect ratio 
changed to 2.0, the results ([triangle] marks) agree well 
indeed with the computation results. If this be so, it 
would seem well to do all the tests with the 2.0 aspect 
ratio, but a high volume of air becomes necessary when 
the Reynolds number becomes high, and data produc- 
+ ee becomes extrernely poor, so this cannot always be 
e. 


The comparisons so far have had to do with such 
localized phenomena on the model surface as shock 
wave position and pressure distribution. Next let us look 
at the drag values where a higher order of precision is 








Re Transition 



































@ NAL, corrected 21x10° free 
(] Harris 9x10" = fixed 
© McDevitt 12.1x10° free 
— Takanashi 21x10° ~— fully turbulent 
-1.6 
-1.2 
C, 
0' 
4 
8 
1.2 




















02 4 6 8 1.0 
x/e 


NACA0012 M=08 a=0° 


Figure 2. Computed Pressure Distribution and Test 
Pressure Distribution—Comparison I 
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Figure 3. Computed Pressure Distribution and Test 
Pressure Distribution—Comparison I! 





required. In Figure 4 we have plotted both the NAL and 
the Harris test data with the Mach number on the 
horizontal axis and the coefficient of drag on the vertical 
axis. The Harris drag is large because the Reynolds 
number is low, while the drag values from the NAL 
wind-tunnel tests are lower than the values calculated 
with the assumption that there is a transition at the 
leading edge, but agree well with the values computed 
when the transition point is fixed at 5%. The NAL tests 
were conducted with natural transitions. 


Next we take up the issue of the reliability of the results 
given by CFD when lift is involved. In Figure 5 we 
compare the theory and the experiments for Mach = 
0.777 and Alpha = 1.536°. Both agree well overall, 
including agreement in shock wave position and inten- 
sity. (In this test we used a model having an aspect ratio 
of 1.5, and the computations were done with the transi- 
tion point fixed at 5% of the wing chord.) However, upon 
closer examination, we see that there is a discrepancy 
between the two at the root of the shock wave. (We have 
already discussed the cause of the difference in pressure 
immediately forward of the shock wave.) This is because 
the shock wave and the boundary layer interfere, so that 
the boundary layer separates at that point and forms a 
small bubble above the wing surface. The Baldwin- 
Lomax is snid to not function well for separation flows 
such as this, and that turns out to be the truth, but if the 
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Figure 4. Mach Number vs. Drag Coefficient 
Characteristics 





flow again adheres to the wing surface, then the agree- 
ment between theory and experiment again becomes 
good, as is evident from Figure 5. This indicates that, 
even if there is separation, the Baldwin-Lomax model 
produces good flow approximations so long as that 
separation is localized at the root of the shock wave. 
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next step naturally would be to do the comparison 
three-dimensional airfoil. Unfortunately, however, 
the three-dimensional case, we have no available 
sion data which can be trusted or backed up by 
pioeoteng ce aatinay ape moh | oF brn 
pressure distributions but no drag or lift inf 
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if it does provide such information it is often very 
difficult to predict correct drag values due to the influ- 
ence of model-support stings or side-wall effects. 
look only at the pressure distributions, it has already 
been demonstrated that the CFD results agree very well 
with the wind-tunnel tests so long as separation is small, 
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O NAL, corrected, Rew3ix10° 
C) NAL, corrected, Rem40x 10° 
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Figure 6. Drag Coefficient vs. Lift Coefficient Curve 








are in a wide variety of ways. 
can be used, for example, |) to pass isobars over a 
three-dimensional wing, 2) to designate lift distributions 
or moments, 3) to suppress undesirable leading edge 
peaks, 4) to weaken shock waves, 5) to prevent trailing- 
edge separation, 6) to achieve wide layer flow regions, 
etc., etc. Many different procedures have been proposed 
to date for resolving such inverse problems, but they can 
be roughly divided into the following three categories: 1) 
Direct methods in which, in principle, pressure distribu- 
tions (or velocity distributions) are designated on bound- 
aries, and then solutions that satisfy these are obtained 
from potential equations or from the Eulerian equations 
30 as to derive the desired airfoil cross-sections. 2) 





terpreting the NACA0002 test data, and to 
Ryoichi Ito of Taiko Denshi Co., Ltd., for his coopera- 
tion in handling the CFD code. 
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[Text] Abstract: With the evolution of the CFD tech- 
nology, the role of wind tunnel has come to be re- 
evaluated. In this paper, the impact of CFD validation 
with wind tunnel is discussed. That is going to bring the 
needs for high accuracy of test data, and that will finally 
lead to more accurate new test technique, e.g., magnetic 
suspension, free from the mount interference of high RN 
low speed wind tunnel, actually very useful for the low 
risk aircraft development. 

Also, wind tunnel facilities in FHI aerospace division are 
introduced. With CFD technology and test technology, 
FHI has. the capability of wide range class airframe 
development; from low speed small target plane to 
space-vehicle. 


As CFD becomes more and more powerful, we now must 
re-examine the role of wind tunnel testing and perform 
revalidations and re-evaluations. In this article we intro- 
duce the reader to some viewpoints of airframe manu- 
facturers to wind tunnel tests and CFD. We hope this 
will enable the reader to better understand how CFD 
research is progressing, in connection with wind tunnel 
testing. We omit here in detailed discussion of CFD 
development and implementation as that was covered in 
the previous symposium. ' 


1, Problems Surrounding Wind Tunnel Tests and CFD 


The questions most frequently encountered in connec- 
tion with CFD and wind tunnel testing are the following 
two. 


¢ Do CFD and wind tunnel testing agree? 

e Will wind tunnels be replaced by CFD? (That is, if 
CFD gets a lot better, will we cease building wind 
tunnels?) 


The answer to the first question is clear, namely that 
CFD which does not agree at all (i.e. does not coincide in 
physical terms) has no technological value, nor is there a 
method for completely computing the motions of air 
molecules. In other words, the issue is the degree to 
which they agree. In answer to the second question, if 
CFD were “almighty” we would not need wind tunnels, 
but such a state of affairs is still far in the future. 
(Landahl, for example, says that, in terms of the calcu- 
lations for reproducing turbulence movements in their 
entirety, where R,, = 10° it would take a year or so, and 








standard problems, but 
when the continuity with the design is considered, 
one ends up carrying on verification checks with wind 
tunnels that are nearby. We cannot have the answers 
agreeing with only one wind tunnel, so we wind up 
doing verification checks using data from a number of 
wind tunnels. This, of necessity, involves us in veri- 
fying one wind tunnel another, and, further- 


wind 
cision is improved during the course of CFD 
cation, in which sense we can hope 
based combination of CFD and wind tunnel 
will result in advances in aerodynamic design tech- 


3. Wind Tunnel Test Reliability 


In Table 1 we represent categorical reasons for attrib- 
uting reliability to wind tunnel test data. We think 
that the concept of “seeming to be reliable” 
comes from the “fact that many tests are conducted, 
which is listed last. We get a sense of all being well 
that does not necessarily agree with the physical 
exactitude because of a number of things, as that 
we are somehow or other able to design aircraft after 
conducting many tests, that we can appreciate the 
continuities with other wind tunnels and understand 
the comparisons with various and sundry derived 
values, and that we obtain data which can be physi- 
cally interpreted, etc. However, when it comes to 
comparison/verification checks on CFD data, phys- 
ical exactitude is in fact demanded, and we come face 
to face with the problems which wind tunnels natu- 
rally have. Typical problems are illustrated in Figure 
1, for none of which have we found determinative 
solutions with Japanese wind tunnels. We cannot 
disregard these problems when we come to consider 
CFD verification checking. We will now discuss each 
one briefly. 
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Table 1. Reliability of Wind Tunnel Test Data 





























Verification of air flow character- | Mach number, sir velocity verifi- 
istics cation 
Static pressure, dynamic pressure 
distribution verification 
Turbulence factor verification 
Drift verification 
Verification by comparison AGARD model 
testing /__ONERA model 
Other standard modeis 
Fact that many tests are conducted 














Figure 1. Wind Tunnel Problems 





i) Wind Tunnel Wall Interference 


change in angle of attack Aa is directly proportional to 
C, we obtain the compensation value by comparing the 
change in pressure integral drag to C, with the change in 
the wake survey values. If the ratio between chord and 
wind tunnel wall height is made small, the interference 
amount can be reduced, but we need a correction in the 

0.6° with C, = 1 even 
is 1:10, which is not a quantity which can 
the effect of the side-wall 
to 


Whenever wind tunnel tests are planned, one must 


pressure (base drag) is very different depending on 
whether a top support or back support is used, it is 
important to conduct final free-flight verifications. In 
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Figure 2. Wall Interference in Two-Dimensional Wind-Tunnel Tests 





(a) Upper Fuselage Sting Support 

















Figure 3. Support Interference Effects on Base Support 





research on aft fuselage flows, the method of supporting 
both main wingtips is sometimes employed (Figure 4), 
that does not mean that it is possible to eliminate the 
interference which develops between the main wing and 
the walls at both tips. Expectations are high for magnetic 


suspension technology in wind tunnels, but it will no 
doubt require computations with models which include 
supporting members in conducting the comparisons of 
exacting wind tunnel tests and computations before that 
technology becomes widely used. 
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ne Results 








Figure 4. Measuring Aft Flows, Example of Results 





The RN effect, as is well known, is particularly pro- 
nounced in flows which include subsonic C,,,,,. separa- 
tion, etc., and the impact on airline design is great. In 
ordinary wind tunnels, however, due to inadequate RN, 
it is not possible to determine the stall characteristics 
and C; max Well, which brings an element of risk into the 
flight tests. In order to verify the effectiveness of CFD in 
aircraft design, it will probably be necessary to verify the 
continuity between actual flight level RN and ordinary 
wind tunnel level RN, but we have no wind tunnels in 
Japan with such high RN. From the perspective of the 
manufacturer, who must fly the actual aircraft, innova- 
tive experimental techniques are employed to make up 
for this RN insufficiency. Here at FHI, since the devel- 
opment of the T-1, we employ a turbulence grid which 
artificially increases turbulence and simulates separation 
conditions at high RN (cf. Figure 5). As shown in Figure 
6, in the case of a straight wing aircraft, the lateral 
instability at stall demonstrated without the turbulence 
grid disappears when the turbulence grid is attached, so 
that we can predict low risk during flight tests. However, 
with a swept-wing craft, the lateral instability at stall 
does not disappear even when the turbulence grid is 
used, and we find out that some kind of device like a 
fence is necessary. Nevertheless, this kind of technique is 
an expedient, and in the future we will need to have 
high-RN low-speed wind tunnels for higher-precision 
development. Also, we do not believe that CFD can 
adequately function in this area at present, and we must 
await further research. 


We can summarize much of what was said above on 
CFD verification as in Table 2. The strict comparison of 
two-dimensional chords and two-dimensional wind 
tunnel tests contains nonsensical elements somewhere, 
but “two dimension” is useful taking design concepts 
from two dimensions to three dimensions, and we 
believe the verification is meaningful in terms of 
grasping the correlation between 2D CFD and wind 
tunnel tests and determining the scope of effectiveness. 





Table 2. Verifying CFD With Wind Tunnel Test Data 
Items To Consider for High-Precision Verification 


2D: Making mode's which include 
wail interference* 


Making 3D models of separation 
eddies* 














RN effect 
Model deformation 


3D: Making models which include 
wall and support interference 


RN effect 
Mc el deformation 
* These two items require three-dimensional models 























4, Wind Tunnel Tests at FHI 


The current status of CFD was treated in the previous 
symposium,’ but we now discuss briefly the wind tunnel 
testing capabilities of Fuji Heavy Industries. We hope 
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Figure 5. High-RN Simulation Using Turbulence Grid 























Figure 6. RN Effect and Leading-Edge Device 
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that this, and the previous treatment, will give a better 

understanding of “wind tunnels and CFD at FHI.” 170 —* 1200 
FHI currently operates six types of wind tunnel and 
water tank with which we can conduct | ao eh AS 
various kinds of tests from low speeds to Mach 4. We = 
believe that we need, at minimum, the equipment to — — 
allow us to perfect basic airframe designs in house RE OER 
——— acerca 

i) With the high-speed wind tunnel (Figure 7 [not repro- | inal 








pore Rete ce -engeiaees dw muy Pmt de pone 
measurement section measuring 0.61 meters square. 
Full-craft six-component tests and various other tests are Table 3. FHI Wind Tunnels 
performed (cf. Figure 4). Since this is a suction type of 
pn gee Rega lad ogg Bang or ag amg 1 pad 
simple, operating rate is high. 11- 
we have conducted some 86 blows with this maximum speed of approximately 80 m/s. This 
raity facility is used in full-craft shape force and pressure 
tests, flutter tests, object dropping tests, ground effect 
ii) The low-speed wind tunnel (Figure 8) measures 2 _ tests, and rotor tests. This wind tunnel is very com- 
m x 2 m, has a normal speed of 50 m/s and a pact for its capabilities. 

















Main Specifications of Low-Speed Wind Tunnel 


Type: Horizontal single-passage circulating type 
Power: 800 KW, 1066 HP DC motor 

Blower: 12-blade, 3.66 m ¢ [diameter] aluminum 

alloy propeller, maximum 650 rpm 
Control method: Stationary Ward leonard type 
Measurement section: 2 m x 2 m, open type 
Wind velocity: Normal 60 m/sec, maximum 80 
m/sec 

Model support: Sting support and strut support 
Data processing: HP-9000 on-line processing 














Figure 8. Exterior View of Low-Speed Wind Tunnel 
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Type: Suction type, two-dimensional wind tunnel 
Blower: 40 HP motor, 8-blade 1000 ¢ [diameter] propeller 
Measurement section: 2m x 2.5 mx 0.2 m 

intake port: 2mx3m 
Wind velocity: Maximum 40 m/sec; 10 m/sec for flow- 








Figure 9. Exterior View of Smoke Wind Tunnel 





iii) The smoke wind tunnel (Figure 9) has a 0.2 m x 2m 
measurement section and uses smoke for flow-line obser- 
vation. It also features three-component, pressure mea- 
surements at wind speeds up to 40 m/s, and is useful in 
the development of airfoils for use at low RN.ꝰ 


iv) The simplified wind tunnel (Figure 10) has a mea- 
surement section that is 0.6 m x 0.6 m and is used in 
visual 3D tests at wind speeds up to 15 m/s. It is useful 
ip tests in which lasers are used to make eddies visible. 


v) The shallow-bottom water tank (Figure 11) simulates 
flows of up to M = 5 in a thin water tank. It is useful in 


quantitative pressure measurements as well as in intu- 
itively understanding the behavior of supersonic and 
hypersonic intakes.‘ 


vi) The visualizing water tank (Figure 12) has a 0.4 m x 
0.6 m measurement section and uses water flows of 50 
cm to 50 cm/s [sic]. It is useful in making the 3D eddy 
sites around the object visible. 


Thus we have various kinds of wind tunnel facilities. As 
for a large-scale wind tunnel such as is needed for highly 
precise development work, however, we must look to the 
installation of facilities such as the aforementioned high- 
RN low-speed wind tunnel in public institutions. 

















Main Specifications of Simplified Wind Tunnel 


Type: Continuous blow-in type wind tunnel 

Measurement section: 60 cm x 60 cm 

Power: 5 HP motor, 4-biade 960 ¢ diameter} 
propeller 

Wind velocity: Maximum 15 m/sec 

_Figure 10. Exterior View of Simplified Wind Tunnel 
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Figure 11. Exterior View of Shallow-Bottom Water Tank 




















x 600 mm, length 730 mm 
Flow velocity: 0 to 50 cm/sec 
Water passage time: Approx. 11 minutes at V = 
5 cm/sec; approx. 20 
seconds at V = 50 cm/sec 
Water passage interval: 30 minutes 
Tracer: Dye or hydrogen bubbles 
Model support:_ Sting support (variable angle) 








Figure 12. Exterior View of Visualizing Water Tank 





5. Afterward 

In writing this article we obtained permission to cite in 
part the results of commissioned contract research, and 
wish to express our gratitude to the Japan Aviation 
Development Association and to the National Space 
Development Agency. 
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[Text] Abstract: At Kawasaki Gifu, the wind tunnel tests 
and CFD technologies have been mostly developed since 
the beginning of 1970. 3D transonic wind tunnel test 
with 1 m x 1 m test section and 3D Navier-Stokes 
Analysis are now available in the design work for many 
airplanes. For example, some airplanes (STOL, T-4, etc.) 
and helicopters (BK117, etc.) were developed by using 
these useful technologies. 


This report presents the current status of these technol- 
ogies in Kawasaki Gifu and the future roles respectively. 








1. Introduction 


Serious wind tunnel testing and CFD work related to 

aerospace got underway at Kewasaki Heavy Industries 

Gifu at about the same time, namely in the late 1960's. 

Currently this work has progressed to a three- 

dimensional transonic wind tunnel and three- 
Navier-Stokes CFD 


During this time, we have made wind tunnel tests and 
CFD the two pillars of aerodynamic design, and have 
employed the two methods in various kinds of R&D and 
aircraft t (modification and new designs), 

with experimental STOL aircraft, the T-4 jet 
trainer, and the BK117 helicopter. We obtained the 
initially performance, and confirmed that the 
technology is effective. 


In this article we discuss the wind tunnel and CFD 
ee 
ture. 


2. Kawasaki Wind Tunnel Tests 


The history of wind tunnel activities at the Gifu facility 
are diagrammed in Figure |. 

Serious wind tunnel testing got underway in 1968 with 
the rebuilding and remodeling of the low-speed wind 
tunnel originally built in 1938. This wind tunnel facili- 
tates both open and closed tests, with the closed config- 
uration being usually employed in helicopter tests for 
rotor blow-down effect tests and the like. Subsequently, 
the two-dimensional smoke wind tunnel was brought on 
line in 1969 and the two-dimensional transonic wind 
tunnel in 1976. The smoke wind tunnel can readily 
capture flow phenomena, including separation flows, 
and can now handle force measurements also with the 
incorporation of a three-component balance scale. The 
two-dimensional transonic wind tunnel is used primarily 
in studies of airfoil aerodynamic characteristics, and was 
put to good use in designing the T-4 airfoils. The 
three-dimensional transonic wind tunnel was built in 
1988 and can be used in a wide range of wind tunnel tests 
from low speed to transonic speed. This wind tunnel has 
a | m x 1 m measurement section, and was designed to 
share in common the models, stings, and measurement 
equipment with the supersonic wind tunnel of the same 
size at NAL. By using this wind tunnel in conjunction 
with the NAL supersonic wind tunnel, tests can be 
conducted continuously over a broad velocity spread, 
from low to high Mach The manifold pressure 
can be raised to 5 kg/cm’, so, for a three-dimensional 
wind tunnel, this provides the highest Re numbers 
obtainable in Japan. The wind-tunnel features are sum- 
marized in Table 1. 


Table 1. Kawasaki Gifu Wind Tunnels 


—Low-Speed Wind Tunnel: Built 1938, rebuilt and 
remodeled 1968 (Open type); Measurement section: 
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“\qure 1. Kawasaki Gifu Wind Tunnel Development 





2.5 m x 2.5 m; Maximum velocity: 60 m/s (Closed 
type); Measurement section: 3.5 m x 3.5 m; Maximum 
velocity: 30 m/s 


—2D Smoke Wind Tunnel: Built 1969; Measurement 
section: 0.1 m x 1.0 m; Maximum velocity: 30 m/s 


—2D Transonic Wind Tunnel: Built 1976; Measurement 
Section: 0.1 m x 0.4 m; Mach number: 0.5 x [sic] 
1.2—2D [sic] Transonic Wind Tunnel: Built 1988; 
Measurement section: 1.0 m x 1.0 m; * number: 
0.2 to 1.4; Max. confluence pressure: 5 atmospheres; 

Re number: 2x 10’; (H = 0.8, C = 0.3M) 


During this period, when broader-region wind tunnel 
en ae ieee ane eee 
wind tunnels and foreign wind tunnels (RAE, NAE, 
and worked hard to raise the level of our 


3. CFD at Kawasaki Gifu 


Kawasaki develops its basic aerospace software at the 

Gifu factory, and from that base also develops applica- 

tion software (for thermal flow analysis and blade/wing 

ee een ane > SE eee oe ee 
visions. 


When CFD is used in computer-based aerodynamic 
analysis and design, we can divide it into the sequential 
problem analysis in which flows around a shape are 
analyzed, and inverse problem analysis in which a design 
condition such as a pressure distribution is given and the 
shape is then derived. We will discuss these two 
approaches while focusing on developments and appli- 
cations in the aerospace field. 


3.1 Sequential Problem Analysis 


Sequential problem analysis came through lift line theory 
analysis and began with the two-dimensional panel anal- 
ysis in 1970, about the same time that we at Kawasaki 
rebuilt the low-speed wind tunnel. It subsequently passed 
through the development of total-craft panel analysis, 
three-dimensional full potential flow analysis, total-craft 
Eulerian analysis, and full-total-craft Navier-Stokes anal- 
ysis to the present stage of obtaining equilibrium flow’ 
analysis for ultrasonic flows (cf. Figure 2). 
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Figure 2. Development of Sequential Problem Analysis Figure 4. Example of Jet Flap Lift Area Theory Analysis 





As to the individual analyses, precision verification 
checks are performed, ing on the correspon- 
dences and so forth with the theoretical values and 
wind-tunnel test results, even after development, and 
improvements are still being made in expanding analysis 
software features, enhancing precision, and I/O data 
processing. 


We now noted some representative analysis results. We 
omit here any treatment of equilibrium flows because 
this is given detailed treatment in Reference 9. 


The panel method is based on non-viscous theory, but 
viscous effects can be considered too with the addition of 
a boundary layer. hid aay Aho ag hy sc 
tion flows by releasing eddies that satisfy the 
conditions for each time step as diagrammed in Figure 3. 
In Figure 4 we have the development to jet flap analysis 
with the addition of aft jet flows, which is employed in 
the aerodynamic analysis of STOL aircraft. The panel 
method, while based on non-viscous theory, exhibits 
good correspondences with wind tunnel tests in regions 
where there is little separation. A , although 
the precision is less than with the differential method, 
the panel method requires less time for computation and 
input data preparation, and we believe it is still an 
effective design tool in parametric studies. ' 


ee) cree a oe ae eS 
comparisons between wind tunnel test results and the 





: # 





results of Navier-Stokes (NS) analysis based on the 
differential and limited volume methods. With these 
methods, it is possible to conduct analyses over a wide 
range, from low to high Mach numbers, including vis- 
cosity. In Figure 5 we compare the calculations and tests 
on the ONERA M6 wing for the purpose of verifying the 
precision of three-dimensional NS analysis. Good agree- 
ment is exhibited, including the position of the shock 
wave. 


In Figure 6 we show an example of the vent effect of a 
ramp and throat in a supersonic intake. An example of 
ground effect analysis is depicted in Figure 7 for a 
vehicle that incorporates a guideway. There is correspon- 
dence with the wind tunnel tests in the upper part of the 
diagram (with the guideway and ground moving together 
with the vehicle), while the lower part of the diagram 
corresponds with an actual vehicle (with the guideway 
and ground fixed, and only the vehicle moving). We see 
differences in separation and force between these two, 
and thus find that care must be exercised in handling the 
wind tunnel test data when inferring actual vehicle 
characteristics from the wind tunnel test results. 


We recently have been conducting studies on STOL 
craft, the HOPE, and the SST, as well as elemental 
analytical studies on flap gap effects.’ 


Thus we see that, thanks to the development of software 
for the total-craft Euler method and total-craft NS 


© Arabyem ~ welts 
<4 Wind-Tunnel Test Results Re = 1.4x10° 
-&: Wind-Tunnel Test Results Re = 2.5x10° 





os 18 4 20 ro 


| = eplitter plate length 
h = bese height 


Figure 3. Example of Eddy Release Analysis 
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Figure 5. ONERA M6 NS Analysis 

















Figure 6. Supersonic Intake NS Analysis Examples 





est Correspondence 

(Guideway Fixed, 
Only Vehicle Moving) 

Figure 7. Vehicle Ground Effect NS Analysis Example 





method, the analytical capabilities have advanced far 
beyond those of the panel method. For this reason, we 
are now able to conduct fine flow studies and do simu- 
lations in modes that are not practicable with wind 
tunnels, although there is still the problem of precision 
when large separation is involved. 


Even with the advances in computers, much time is 
reauired from innut data prenaration to the finished 


and these problems need to be resolved before this 
method can be used in the design stage. 


3.2 Inverse Problem Analysis 


In design work in which wind tunnel tests and CFD 
sequential problem solving is used, shape modifications 
and aerodynamic analyses are often repeated. With this 
approach, however, there is no guarantee that the design 
goals will be reached, and a great deal of time is required 
to obtain convergent solutions. 














that of the design of a main wing which includes the USB 
engine. 


4, Wind-Tunnel Tests and CFD 


In Table 2 we compare wind-tunnel tests and CFD in 
terms of various problematic categories. 


Table 2. Wind-Tunnel Tests (WTT) and CFD 
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—Reynolds number, Mach number: WTT: Limited 
(especially Reynolds number); CFD: Precision is prob- 
lematic, but yields solutions that correspond to Re 


numbers, Mach numbers of actual craft 

—Flow Simulation (separation flows): WTT: Can obtain 
actual flows (yields solutions that include stall); CFD: 
Limited with respect to physical model com- 
prehending large separation is very | ) 

—Interference Problems (air flow x WTT: 
Wind-tunnel tests are very difficult; Can obtain 
ideal solutions 


— e eg hgh ted 


differential method or limited X method 
—Reproducibility: WTT: Depends on model precision, 
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Figure 10. Three-Dimensional NS Inverse Problem Analysis Example 





with sophisticated design technology and instruments. 
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understand the degrees of approximation and simulation 


which CFD provides, so as to appreciate the scope of 
applicability. 
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5. Concluding Remarks 


The biggest problem with wind tunnel tests is the failure 
to achieve actual aircraft Re numbers with wind tunnel 
Re numbers. For this reason, in order to grasp actual 
craft characteristics, one must rely on Re number extrap- 
olations, with which accurate estimates are very difficult. 
This being so, as can be seen in Boeing’s keen interest in 


On the other hand, there are still problems with CFD 
sequential problem analysis, such as the ability to ana- 
lyze flows which include large separation, the long TAT 
times, and solution reproducibility. We now need to 
re-evaluate turbulence models, shorten preprocessing 
times using non-structural grid techniques, and study 
dependencies using grids. 


We are working in the area of non-structural grid tech- 
niques, having perfected software for two-dimensional 
NS analysis and working on the application of this to 
three-dimensional NS analysis. 
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Snape: Abstract: Computational Fluid Dynamics (CFD) 

has established its unique position as a design tool in 
aerospace industries. As the CFD technology has 
evolved, the relationship with the conventional design 
tool of wind tunnel testing has also changed accordingly. 
This paper reviews the relationship of the two major 
technologies in Mitsubishi Nagoya. The typical CFD 
application is briefly introduced during the discussion. 


1. Introduction 


The utility of CFD technology in the aerodynamic design 
of aerospace equipment began to make rapid strides in 
the 1980s. Mitsubishi Heavy Industries (MHI), at its 
Nagoya Aerospace Systems (NAS) facility, has also been 
working hard on CFD R&D and now considers it an 
important aerodynamic design tool alongside of wind 
tunnel testing. In this article we trace the development of 
CFD at Mitsubishi NAS, giving specific examples and 
focusing on how its relationship with wind tunnel testing 
has changed. We also discuss the significance of these 
two technologies to engaged in aerodynamic 
ae ten tae Sa opal ada 
no 


2. Changes in CFD and Wind Tunnel Testing 


Advances in CFD applications in the aerospace field 
were tremendous in the 1980s, and CFD has now devel- 
oped into an indispensable aerodynamic design tool that 
ranks alongside conventional wind tunnel testing. As 
CFD technology has naturally its relationship 
to wind tunnel testing has changed (cf. Figure 1). 


In the early 1980s, when we began to see the possibilities 
for CFD applications, it was being said that CFD was a 
“substitute” for wind tunnel testing. At that point in 
time, that being focused on in CFD technology were the 
analytical functions, as was the case with wind tunnel 
tests, which put CFD in a competitive relationship with 
wind tunnel testing. Subsequently, with the unforeseen 
improvements in supercomputer speeds and CFD tech- 
nological advances, it became possible to handle Navier- 
Stokes inverse solution computations, and after the 
mid-1980s it was the design functions of CFD which 
came into focus. These were absolutely new functions 
which are not available in conventional wind tunnel 
testing, and in that sense CFD technology transcended 
wind tunnel testing. Then in the late 1980s the quality of 
CFD technology began to be questioned, which led to 
careful comparative studies between CFD and wind 
tunnel tests. As a result, there is now emphasis on the 
complementary functions realized by organically com- 
bining the two technologies. 


Beginning in the next section, we introduce the reader to 
some typical examples of CFD applications at Mitsub- 
ishi NAS, and retrace the changes in the relationship 
between CFD and wind tunnel testing. 
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Figure 1. Changes in Relationship Between CFD Tech- 
nology and Wind Tunnel Testing 





3. Competition (Analytical Functions) 


In terms of aerodynamic analytical functions for deter- 
mining aerodynamic characteristics for a given shape, 
there is no essential difference between wind tunnel 
testing and CFD. However, in the field of aerodynamic 
design, these two technologies are clearly differentiated 
in their use. CFD is superior to wind tunnel testing in 
terms of its greater speed and its greater flexibility. 


Compared to the task of designing and fabricating 
models and then conducting wind tunnel tests, the time 
involved in generating CFD meshes is extremely short. 
For this reason, CFD has great merit in being able to 
immediately reflect the results of aerodynamic perfor- 
mance evaluations in shape designing. And unlike wind 
tunnel tests, where the test conditions are limited by the 
capabilities of ground-based equipment, CFD tech- 
nology in principle poses no restrictions on the evalua- 
tion conditions, so that aerodynamic performance eval- 
uations can be conducted in regions that cannot be 
simulated with a ground-based facility. In Figure 2 a 
pressure distribution is depicted for HOPE re-entry 
conditions. This kind of function is particularly appre- 
ciated in the development of such space shuttle 
vehicles.’ In Figure 3 [not reproduced] is displayed a 
pressure distribution around the blades of an advanced 
turbopropeller (ATP). Obtaining this kind of detailed 
data for a turning object from wind tunnel tests is 
extremely difficult. 


In the area of “analytical functions,” where it is essen- 
tially no different froin wind tunnel testing, there are still 
many respects in which CFD takes a back seat to wind 
tunnel testing. Take the cost angle, for instance. For 
some shapes, where many data points are required, wind 
tunnel testing holds the advantage. And in the area of 
precision, beginning with turbulence models, some 
aspects of CFD are still unknown. And there is a tremen- 
dous difference in terms of track record. In the field of 
aerospace equipment development, for example, where 
very high reliability is demanded, nothing has ever flown 
yet with only CFD. 








Figure 2. HOPE Pressure Distribution (M = 15, a - 
30°, Re = 2.8 x 10°) 


JPRS-JST-94-002 
7 February 1994 





re 
ry 


Figure 4. Full-Craft Pressure Distribution for Transonic 
Civilian Aircraft (M = 0.82, a = 1.5°, Re = 2 x 10°) 





4. Superiority (Design Functions) 


CFD design functions perform processes which are 
exactly the reverse of the analytical functions, and here 
we can expect little help from wind tunnel testing. The 
ability to derive a geometrical shape given a pressure 
field, in what is called inverse solving, is extremely 
important to the aerodynamic designer. Recently it 
became possible to perform Navier-Stokes reverse solu- 
tions, and the field of applications has spread. To the 
aerodynamic designer who must design “flows,” it is a 
tremendous advantage to have the “directness” of being 
able to designate not shapes, but pressure distributions. 


In Figure 4 is depicted a full-craft pressure distribution 
for a transonic civilian aircraft.? The high-Reynolds 
number main wing that is employed in the latest air- 
frames represents the positive application of the Rey- 
nolds number effect—which used to be difficult to 
estimate in many cases—so as to achieve high drag- 
divergence Mach numbers. The technology for these new 
wings would not be possible without the high-precision 
reverse solution method of CFD technology. In Figure 5 
is depicted an example of applying the Navier-Stokes 
reverse solution method.* The HOPE main wing has a 
complex flow field that is surrounded by the fuselage and 
chip fins, so its design required a tool for precisely 
estimating the positions and conditions of the shock 
wave and separation, particularly at transonic speed. By 
employing the Navier-Stokes reverse solution method, it 
became possible to make extremely high-precision 
designs in a short time. 


5. Complementing (Organic Merging) 


It is impossible [sic] to verify the precision of CFD 
computations by doing careful correspondence studies 
with wind tunnel test data. Out of the process of carefully 
comparing the two kinds of simulation results, we derive 
“organic merging” with which it is possible to make 








Figure 5. HOPE Main Wing Reverse Solution Design 
(M = 0.90, a = 5°, Re = 2 x 106) 





aerodynamic performance evaluations at higher dimen- 
sions by mutually making up for each method’s short- 
comings. 


In Figure 6 we have a pressure distribution for a swept- 
forward wing. With a swept-forward wing, due to the 
structural divergence characteristics, the deformations 
that result from static aeroelasticity, even in a wind 
tunnel model, cannot be disregarded.° For this reason, 
we combine CFD and structural analysis to estimate the 
flexure in the model beforehand, and then actually 
fabricate the model in a deformed condition so that it 
forms the aerodynamic shape designed at the design 








JPRS-JST-94-002 
7 February 1994 





Figure 6. Forward Swept Wing Pressure Distribution 
(M = 0.80, a = 1.0°, Re = 2 x 10°) 
































points in the wind tunnel tests. This enables us to verify 
aerodynamic characteristics in a way that 
could not be done with either CFD or wind tunnel 
testing alone. In Figure 7 [not reproduced] are depicted 
the results of two-dimensional wind tunnel side-wall 
analysis.° In two-dimensional wind tunnels, particularly 
in developing transonic airfoils with high Reynolds num- 
bers, it is not possible to use models with high aspect 
ratios due to limitations in model structural strength, so 
the side-wall effect becomes too large to be disregarded. 
For this reason, we employ CFD to analyze the side-wall 
effect in a two-dimensional wind tunnel, then apply 
corrections to the wind tunnel test results, and thus 
perform high-precision performance evaluations on 
transonic airfoils with high Reynolds numbers. 


Now that CFD technology and wind tunnel testing have 
reached the mutually complementary stage, from the 
aerodynamic designer’s perspective the situation 
becomes as diagrammed in Figure 8. Wind tunnel tests 
provide CFD technology with a steady stream of verifi- 
cation data for enhanced precision, while CFD provides 
wind tunnel testing with images of flow-site phenomena 
inside the wind tunnel which is very difficult to see in a 
conventional wind tunnel test. The designer continually 
learns physics from the wind tunnel tests while learning 
how to narrow his focus on the points which should be 
looked at in terms of designing and testing. 


6. Summary 


Amazing progress is being made in CFD technology at 


Mitsubishi Heavy Industries’ Nagoya Aerospace Sys- 
tems facility. In conjunction with these advances, the 
relationship between CFD and wind tunnel testing is 
changing. Aerodynamic designers should now be able to 
effectively use these two tools in a mutually complemen- 
tary way to design higher-performance aerodynamic 
shapes. 
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Figure 8. Two Design Tools 
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[Text] Abstract: When aero-engines are developed, sev- 
eral different types of wind tunnel tests are involved. For 
instance, one type of tunnel is required for fundamental 
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close replicas. Certain effects may not be adequately 
addressed in wind tunnel tests, for example, the effects of 
high speed revolution, flow properties in unmeasureable 


of numerical wind tunnels parallel those of experimental 
wind tunnels. The paper also assesses future trends for 
the use of numerical wind tunnels. 


1. Introduction 


Wind tunnel testing in aircraft engines involves widely 
diverse configurational elements, wind tunnel scale, and 
computation volumes, from the basic turbine and com- 
pressor blade cascade tests to the high-altitude perfor- 
mance tests of the complete engine in the wind tunnel. 
Usually wind tunnel tests and CFD are used together in 
blade cascade tests and duct tests, but the scope has 
broadened in recent years with the phenomenal 
advances made in computers and analytical techniques. 
We here look at the different wind tunnel tests used on 
aircraft engines and discuss how CFD fits into the 
picture. 


2. Wind Tunnel Testing of Aircraft Engines 


2.1 Types of Wind Tunnel Tests 


As diagrammed in Figure 2-1, the wind tunnel tests used 
on aircraft engines include high-altitude performance 
tests (Figure 2-2) to determine overall engine perfor- 
mance, ground performance tests (Figure 2-3 [photo not 


reproduced]), element tests on fans, compressors, com- 
bustors, and turbines (Figure 2-4), and configurational 
element tests such as blade cascade tests and strut tests. 


Then these configurational element problems are further 
divided for investigation into wingtip clearance tests, 
cooling hole effect tests, and the like, and we also do 
basic testing on aft surface step flows and shear flows, 
etc., in analyzing basic internal flows. 


The ultimate purpose of wind tunnel testing is to 
improve the performance of the whole engine. For this 
reason it is not enough for these tests merely to be 
reliable. There must also be a structural rationale to the 
entire series of tests, and they must also match the 
development resources and time allowed. 


2.2 Correspondence With CFD 


Compared to the wind tunnel test facility, which ranges 
widely in size, the test sample [i.e. model] is almost 
always made of the same material and in the same shape 
to withstand high speeds, high temperatures, and high 
pressures, Because this is so, it is more important that 
there be an equivalency between the actual craft’s oper- 
ating conditions and the wind tunnel tests than that of 
the actual craft and the test sample be similar in the 
flight tests and wind tunnel tests. The issue of equiva- 
lency in wind tunnel tests includes such things as simi- 
larity in the supersonic and hypersonic wind tunnel tests, 
and in the environmental impact tests for noise and NO, 
generation and propagation. However, if we limit the 
problems to the flows inside the engine, then the issue of 
the degree to which the wind tunnel inlet/outlet and side 
walls coincide with the actual operating conditions 
becomes very important. Even when we break the whole 
engine down into element tests, blade cascade tests, and 
basic tests, setting the test conditions is still a very 
important issue. When we consider the correspondence 
with CFD, we try to achieve correspondence primarily in 
element tests and blade cascade tests, depending on the 
scale of the computer and test equipment, and the time 
allowed for the wind tunnel tests, and we are working 
very hard to improve the precision with which we set the 
boundary conditions. 





JPRS-JST-94-002 47 



















































































7 February 1994 
Readily Removable 
Cooler | | 
Stage 3 > 
— 
Inlet Control 
Valve 
Royal Attitude Engine Test Facilities Comparable 
Alroraft Facilities 
Establishment 
Oyestock Test Chamber Size: Mase Flow: DOD-AEDC. 
nites (ft) 25 dia x 56L | (b/sec) 1400 | +1, 42,0-1,C2 
Kingdom Paw: %-217, X-218 
Date Bullt/Upgraded: Altitude range: 
1969 | () 50,000 
ATF Replacement Cost: Temperature Range: 
Cell 3W £60M* (F) -60 to ambient 
Operational Status: Pressure 
Double day shift (pela) 2 to atmospheric 
Speed range 
(Mach No.) Subsonic 
Direct-connected and freejet, high-acouracy thrust Group 4 
capability: 6000 Ib/1 
Special cepability le connected or free loing. 

















Figure 2-2. High-Altitude Performance Test Facilities 
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Figure 2-4. Fan Element Test Apparatus 





3. Roles of Wind Tunnel Tests and CFD 


1) Prediction of R&D results; prediction of problems, 
possibilities prior to design, testing 


2) New design, modification design; designing by adding 
CFD technology—which analyzes important design 
points with suitable means in the limited time allotted— 
to conventional means 


3) Numerical/mathematical wind tunnel testing; nar- 
rowing the focus of wind tunnel tests and substituting for 
wind tunnel tests 


The object of analysis will vary with the size of the 
computer, the reliability of the analysis, the tire 
allotted, and the costs allowable; and we think it will take 
on the same type structure as the wind tunnel tests 
described above. 


3.2 Predicting RAD Results 


In Figure 3-2 we have diagrammed the prediction of 
R&D problems and possibilities by wind tunnel testing 
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Figure 3-1. Role of CFD 




















Figure 3-2. Prediction of Problems, Possibilities 


and CFD. With correspondence to wind tunnel tests, 

CFD is being used to verify the presence or absence of 

separation in flows generated by simple blade cascades, 

and in elucidating pronounced well interfere 
num 


gradients as in fans and compressor blade cascades. In 
Figure 3-9 we show a comparison in the case where the 
injection angle is varied, inclusive of non-design points. 


~ —* 
ey 4 
=> ce 


— a 
fm ie 


% 

















Figure 3-5. Flow Velocity Vectors Inside Combustor 





We believe that this has become useful in estimating 
qualitative trends, particularly injection angle limita- 
tions. In Figure 3-10 we give an example where these 
were used in a modification design, and in Figure 3-11 a 
comparison with wind tunnel test results. By these 
means we can verify the propriety of both our predic- 
tions and our modification designs. 


What is urgently needed now is the development of 
techniques for analyzing blade cascades which involve 
secondary flows in reverse pressure gradients. 


3.4 Numerical Wind Tunnel Tests 


We now need numerical wind tunnel tests which will be 
nearly as reliable as, and serve in part as substitutes for, 
conventional wind tunnel tests. To create such numer- 
ical wind tunnel tests, we must perform more of the 
verification checks discussed earlier, do detailed studies 
on turbulence models and schemes, and improve our 
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Figure 3-6. Comparison of L2F Measurement Results and Numerical Analysis Results 





























Figure 3-7. Comparison of L2F Measurement Results and Numerical Analysis Results 
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Results of 3D Flow Analysis on - 
Static Fan Blades for IHi Research 
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Figure 3-8. Comparison of Static Fon Blade Analysis Results and Experimental Results 
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Figure 3-9. Variation in Flow-Out Angle and Pressure Loss When Varying Injection Angle 
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Figure 3-11. Modified Design Wing Aerodynamic Performance Characteristics Comparison 














flows that cannot be conventionally measured. 
ingly, this foundational measurement of internal flows 
Neen Crt ee ee oe 
estimations of internal flows. In the past the focus has 
been on measuring the pressures on wing surfaces and in 
blade cascade inlets and outlets, but the scope of appli- 
cations is broadening to such things as the photographic 
effects of the behavior of shock waves at the leading 
edges of rotating fans by means of holographic measure- 
ments, depicted in Figure 4-1, and the measurement of 


these will pave the way for improved CFD reliability and 
wider CFD applications. 


5. Concluding Remarks 


The scale of wind tunnel tests in aircraft engines varies 
widely from high-altitude performance testing on the 
high end to basic blade cascade tests on the low end, 
ee a oe OO ere oe 
re ety Gee oa ood p thoy tg) | 


importance than wind tunnel testing and being 
mutually complementary in high Reynolds number 
regions, namely how CFD can be wed to etimate wind 
tunnel tests with good precision, providing inf 
which cannot be measured in a win 

reflecting this in the design. And, as CFD 
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Changes in Shape of Shook Wave Between 
Wings Due to Stopping Down (96% Rotation) 


Figure 4-1. Inter-Wing Shock Wave Shape Based on 
Holographic Measurement 





give numerical wind tunnel 


breakthroughs 
testing a more important role to play. We would like to 
take up the matters of environmental impact survey 


testing and supersonic/hypersonic testing in a subse- 
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